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We study the distributions of the continuous-time quantum walk on a one-dimensional lattice. In 
particular we will consider walks on unbounded lattices, walks with one and two boundaries and 
Dirichlet boundary conditions, and walks with periodic boundary conditions. We will prove that all 
continuous-time quantum walks can be written as a series of Bessel functions of the first kind and 
show how to approximate these series. 
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I. QUANTUM WALKS 

Quantum walks have been introduced as an algorith- 
mic idea for quantum algorithms in an attempt to speed 
up classical random walk algorithms, see e.g. 0, ITT| . 
Two different proposals were made: the discrete-time 
quantum walk [2, [l3| and the continuous-time quantum 
walk [i3,|f|- While the behavior of the discrete-time quan- 
tum walk is well understood by now, the continuous-time 
quantum walk has not been studied as extensively, but 

see Hi mm. 

In this paper we consider the one dimensional con- 
tinuous-time quantum random walk on a line. We write 
the state of the random walk system after time t as 



kh*)> = X>(m) i*) 



where L, R 6 Z U {±oo}, L < R, are the left and right 
boundaries of our walk. Our goal is to determine ip(x,t) 
in an analytic form. A continuous-time quantum random 
walk is now defined by the difference/differential equation 



*9i 



-ip(x-l,t)+qip(x,t)-ip(x + l,t) (1) 



and boundary conditions on ip(L, t), ip{R, t). 

In this paper we study four kind of boundary condi- 
tions: the unbounded walk (L = — oo and R = +oo), 
walks with one boundary (L finite, R = +oo, and 
^){L,t) = 0), walks with two boundaries {L, R finite, 
ij)(L,t) = ip(R,t) = 0), and periodic boundary condi- 
tions (L, R finite, ip(L,t) = ip(R,t)). 

If both L and R are finite we define N := R — L as the 
"length" of the interval in which the walk takes place. In 
this case the continuous-time quantum walk from JQ) can 
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also be written as 



4 t \m) 



-H q m)) 



(2) 



Here H q is a (N — 2) x (N — 2) matrix, when we use the 
ip(L,t) = ip(R,t) = boundary conditions, or a (N — 
1) x (N — 1) matrix for periodic boundary conditions, 
which is defined as 



H„ 



-q 1 
1 -q 1 



-q i 
i -q. 



(3) 



on the vector space span{|L + 1) 
for i/)(L,t) = 4>(R,t) = 0, and 



|£ + 2),...,|i?-l)} 



Hq tV — 



q i 



-q 1 

l -c 



(4) 



\R- 1)} for 



on the vector space span{|i) , \L + 1) , 
periodic boundary conditions. 

The solution | of © for a given starting state 
|V>(0)) is given by 



iH„t 



(5) 



which gives yet another definition of the continuous-time 
quantum random walk. 

In this paper we will prove the following theorem: 

Theorem 1. Let J n (t) denote the n-th Bessel function 
of the first kind, see e.g. Q/. Define J n {z) — i n J n (z). 
The continuous-time quantum random walk defined by 
equation QJ) with starting distribution |?/>(0)) = \xq) has 
the following coefficients ifj{x,t) at time t: 

1. For an unbounded walk, i.e. L = — oo, R = +oo, 



^{x,t) = e-^J\ x _ X0 \{2t). 



(6) 
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2. For a left boundary LeZ and R = +00 and bound- 
ary conditions tp{L,t) = 0, 



*p(x,t) 



J\x-x \(2t) - J\x+x -2L\{2t) 



(7) 



3. For two boundaries L, R € Z, L < R and boundary 
conditions ip(L, t) = ip(R, t) = 0, 



1>(x,t)=e- it * J2 (-i) n J\x-x nl m, 



with 



L — (x- n -i — L) n < 
R + (R - x- n+ i) n>0 



(8) 



(9) 



For two boundaries L, R € Z, L < R with periodic 
boundary conditions ip(L,t) = ip(R,t), 



i/j(x,t) 



- itq 



E J \x-Vn\( 2t )' 



with 



y n := nN + x . 



(10) 



(11) 



The unbounded case ((HJ) was known before for q = 
0,2, see 0,0. For q = equations and JHJ were 
derived by Apolloni and de Falco in their treatment of the 
clock of a quantum computer |4| . The walk with periodic 
boundary conditions ljl0|l was studied approximatively in 
0, E| . We will provide a unified proof of all these results 
in this paper. 

The infinite series (jHJ), i|10[l are of course unwieldy for 
practical applications. We can truncate JHJ) and l)10|) af- 
ter a certain number of terms and still get a very good 
approximation. 

Theorem 2. For two-sided and periodic boundary con- 
ditions the series 0) and l^l(J\) . truncated after 



k = 



t + N 21nt+ M<T 1/2 ) 



N 



ln(21nt+ '" (cE -r l/2) ) 



terms, with constant c :— ^| , approximates ip(x,t) up 
to e, i.e. 



tp(x, t) 



(-l) n J\x-x n \(2t) 



< £ 



, l . i ^ N _ l 

for < e < 1 and t > max{l, I e e (£/c)« ) }. 

We use Theorems ^ El to create a plot of the proba- 
bility densities for a quantum walk. In Figure ^ you see 
probability density plots that were made by truncating 
at the appropriate k given by Theorem [21 



II. ANALYSIS OF WALKS WITH NO 
BOUNDARIES AND ONE BOUNDARY 



We can easily check that for k = 1, . . . , N — 1 

/ feTrf.T:— \ , 

v ^ / k: ' 

x=L+l 

are normalized and orthogonal eigenvectors of the ma- 
trix H q defined in J3J with eigenvalues 2cos(-^) — q. 
Therefore we can write \ip(t)) = e ltHq | -0(0)) as 



N-l 



~L)\ ^it2 cos "jy" J — 



fc=l 



x J2 sin(^pM)^(y,0). 

!/=£+! 



Suppose we start the random walk at xq, i.e. our starting 
state \ip(0)) = \xq) or equivalently i/j(x,0) = S XjXo . Then 



2e~ itq 
N 



N-l 



E sin( " fck L) )sin( 



7rk(xQ — L) \ 
iV / 



fc=l 



x g^cos^ J 



N-l 

E 



A'TT N 



/ 7rk\x — xq\ \ it2cos \ —r- 

cos[ — ' N ' je \n 



k=l 
N-l 



-itq 



N 



^>os( ^ + *°- 2L| y 



iV . 



fc=l 



(12) 

We replace the sums with integrals. The error for these 
approximations are (see e.g. the errors for Newton-Cotes 
quadrature rules in | l0| ; note that the missing endpoints 
k = 0, N can be added without penalties, since they can- 
cel out in O): 



N-l 



fe=i 



f(x)dx 



<0^llr(*)IL 

where the norm indicates the supremum of the second 
derivative of / with respect to x taken in the interval 
[0, 1]. Since for a, (3 > 

2 



dw 2 



cos 



[aw) 



^i(3 cos w 



<(a + (3Y+f3 



we can rewrite equation l|12|) as 



V>(x, t) = e' ltq / cos(7r|x - x \w)e it2cos ^dw 



e~' ltq / cos(Tr\x + x -2L\w)e it2cosi ™Uw 



± e(x, xq, t, N) 
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with an error term of order 



e(x, X0 ,t,N) = o ( a— °i+2*) 2 +a^o-2L| +2t) 2 ) _ (13) 

This representation allows us to express ij){x, t) as a sum 
of Bessel functions. According to pj, equation 9.1.21, for 
n £ N + the Bessel function J n {z) is given by 



J n {z) = / 

I" JO 



JIZ cos w 



cos(nw)dw 



and therefore J n {z) — ^ e lzcosw cos(nw)dw, and we 
can write 



ip(x,t) 



-itq 



u \X — Xq\ 

(2t) - J, (2t) 



± e(x, so, t, N). 



(14) 



For _R — > oo , which implies iV = i? — L — > oo, the error 
term e{x,XQ,t, N) defined in (|13|) converges to and we 
get 



- i I q 



x-\-xq— 2L\ 

(2i) 



which proves J7J|. This is the limiting distribution for a 
continuous-time quantum random walk with one bound- 
ary. For no boundaries, i.e. if we let L — > oo, too, Q 



simplifies to 



1>(x,t) =e- Uq J {x _ Xol (2t), 



since linin^oo J n (x) = 0, which proves (J7J) and can be 
compared to the results from ||J • 



III. WALKS WITH TWO BOUNDARIES 

To derive the correct form of the quantum walk dis- 
tribution for walks with two boundaries we check the 
correctness of Equation © for the unbounded case by 
plugging it into (Q: 



dt 



where we use that the derivative of the Bessel function 
satisfies J' v (z) = J '-' (z) ' J ' +lW , see e.g. Q, to verify 
that ip indeed fulfills QJ, at least for \x — xq\ > 1. In the 
special case \x — xq\ — we use that J' {z) = — Ji{z). 
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To satisfy the boundary conditions il)(L , t) = ip(R, t) = 
we use the method of images, see |l^, and introduce 
mirror images of the walk starting at xq. Let \ip Xo (t)) 
be the solution for an unbounded walk starting at xq. 
Then \ipjkj(t)) = IV'xoW) is a good approximation to 
the solution, but we have to correct the error on the 
boundary by introducing walks starting at x_j = L — 
(xq — L) and x% = R + (R — xq): 

|^ ) (t)> = K ) (*))-|^-i(*))-|^x(*))- 

To correct the error introduced by the corrections we 
continue to add mirror images at x_2 = L — (x± — L) and 
x 2 = R + (R — X-x), until we arrive at iJSJ. Let us 
check that 



rl>(L,t) 



- itq 



E (-l) n ^-»|(2i) 

n— — oo 
oo 

El" 1 )" [j\L- Xn m - J\L- x _ (n+1) \(2t) 

oo 

El- 1 )" \j\L- Xn pt) - J\ Xn _ L \{2t)\ = 0, 



and similarly ip(R,t) = 0. This proves ©■ 

For periodic boundary conditions ip(L,t) = ip(R,t) we 
also build a solution starting from \ipx (t))- I n order to 
enforce periodicity, we add two walks at y_i = x a — N 
and yi = xq + N, which is the correction term for the 
"first round" . If we define 

Vn =x Q +nN 

we arrive at i|l(J[l . Then 

i>(L,t) 



-itq 



= zl J \L-y n \i 2t ) = E J\L-nN-x \{2t) 



71— — OO 
OO 



71 — — OO 



i>(R,t) 



-itq 



= E J\R-{n+l)N-x Q \{2t) 
n— — oo 

which completes the proof of Theorem ^ 



IV. APPROXIMATIONS 

For practical purposes the infinite series in Theorem 
n are unsatisfactory. Theorem [3 which we will prove 
in this section, fortunately tells us that we can truncate 
those series and still get an excellent approximation to 
the correct distribution. 

Consider the distribution ip(x,t) in JSJ for zero bound- 
ary conditions 



ip(x,t) 



OO 



(-l)"J| x _ Xn |(2t), 



and truncate it after the first k terms. We would like to 
determine the error of this approximation, i.e. 



efc.o 



i>(x,t) - e- u " £ (-l) n J|,-,„|(2t) 



n— — k 



for t > 0. We can bound the Bessel functions J u {2t) by 

IC 



l J ^r(, + i)< 

see 0. Then the error e^o is bounded by 

oo — oo 

efe,o S E |^|x-x„|(2*)|+ E | J I 

X — X n \ 

(2t) 

n—k-\-l n——k— 1 

^l^ - x n\ + \x— X—n\ 



(15) 



* E 

n=fc+l 



r(|as - ar„] + 1) r(|x - z_ n | + 1) 



We similarly define the k term approximation error ek }P 
for periodic boundary conditions. 

From the definition 10 of the reflection points x n 

{L — (x_ n _i — L) n < 
R+{R-X- n+1 ) n>0 

we see that L + nN < x n < R + nN. Therefore, and 
since L < x < R, 



oo 

efc,o < E 

n=k+l 

oo 

* E 



£|x-x„| 



£|x-x_„| 



r(|a:-a; n | + l) T{\x - x_„| + 1) 



r((n-l)JV+l) r((n-l)iV + l) 



oo , ni y oo 

<2t w V - - 2t 2N V 

n=k+l vv ; ; n=k 



(rciV)! 



(16) 



The same error estimate also holds for the error ek, p of 
the fc-term approximation to the series (|1L)|) for periodic 
boundary conditions. 

To estimate e^o and ek lP in (|16f) we can use Stirling's 
formula 

x\ = V2^x x+ ie~ x+ ^ > V2^x x+ h~ x , 
see 0, where x > and < 6 < 1. Then 

(niV)! >V2^(nAf)( nJV )+5e-( niV ) 
^V^NinNe' 1 )^. 

Therefore the error e/^o from l|lt)|) is bounded by 



2JV 00 



2t 



V2irN v a 



Y — 



te 
nN 



^N 



< 



2t 



2N r t e ~\ nN 



V2^Nk 



E 



nN 
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This series converges if k = k(t, e) is chosen such that 
a n / te \ Nk °° / te \ N ( n ~ty 



k > |l. In this case 



n—k 



Etc ( I ( 



te 



Nk 



1 



7 1 _ / ^ \ 

7 1 I fcjV / 

-Wfeln 



AT 



(17) 



1 - Ms-V 

Therefore the error from truncating after fc terms is 

efc.o < 



v / 2vre7Vfc l - ( ") 

V /civ / 



A' 



(18) 



We want to use (fT5fl to find k such that efc,o = 0(e). 
We make the following ansatz 

k =k{t,e) = —^1^71 where ( 19 ) 

C=C(*,e) = 21nt+i^i^ 1 (20 ) 

with a constant c > that we will determine. We have 
to check that the requirement k > 4§ that guarantees 
convergence in (|17|l is fulfilled. If we restrict t to 

* > (e e (e/c)Tr 



1 \ 4JV-1 



then 



C = 21nt+ >->e. 



(22) 



Additionally we observe that ^S- is growing monotonous- 
ly for z > e, 



d z In z — 1 



> for z > e, 



dzlnz (lnz) 2 

and therefore for t chosen as in <|21[1 

t + N ( t + N te 
hiC " N~ e > ~N' 



(23) 



We estimate the numerator and denominator of e^o in 
(|T5jl separately and write e^o < wct- If we insert ijT§|) 
into the numerator of (|18fl . we see that 

M{k) : =2t 2N e- Nkla ^ = 2^ e -^^M 1 #^£) 

/ In t + W lnlnC" 

=2^ e - (t+jVK ( 1+ ^ 
Now we use that we can show that for z > 1 
ln^-lnlnz l n *±^-lnlne^ 



lnz l ne ^ t + TV" 

We can use this to estimate M(k) since C > e > 1, see 
(EH, and get 

M(k) < 2t 2N e~< N = 2t 2N e~ 2NXat e x ^ 1 ^ = -eV~t. 



To estimate the denominator of (|18fl 



D(k) := V2neNk [ 1 - ( — 
1 1 fciV 



we insert (|19|l in 
te 

kN = 



te 



t elnC 



LtfS_ N t + N C 



(24) 



We can see that both and 2j2±& are less than 1 for 
t as in l|21|) . when we have C > e, see (|22H . We analyze 
(|24|l by considering t £ [(e e (e/c)~^ ) 4N -! , 00) over two 
intervals. For t S [(e e (e/c)"^)™^i , 10) 



t e In C 



< 



< 



i + iV ~ t+1 



< 



10 

IT' 



In the interval [10, 00) we use that 2J^_k converges 
monotonously to for £ = C(t,s) — > 00, and that 
C > 21n*: 

t elnC elnC eln(21nt) 10 

< < < 0.902 < — . 

t + N C ~ C " 21nt " " 11 



(21) Therefore we can use this and ill' Ml) for 



I e\/2irt 

D(k) > V2TreNk— > . 



If we combine the results for numerator and denomi- 
nator, we can estimate 



M(k) 

ein < — — - < 

'° - D(k) ~ 



22 



ce\r2/K 



and we see that we have to choose c := ? . This result 

e\/27r 

holds for all t chosen according to (|21|l . 

If we follow the same steps we can also estimate the er- 
ror ek, p of a k term approximation for periodic boundary 
conditions. This proves Theorem [21 
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